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In this Letter we discussed the parametric instability of texture of homogeneous (in time) spin precession, 
explaining how spatial inhomogeneity of the texture may change the threshold of the instability in comparison 
with idealized spatial homogeneous case, considered in our JETP Letter 83, 530 (2006). This discussion is 
inspired by critical comment of LA. Fomin (cond-mat/0606760 ) related to the above questions. In addition 
we considered here results of direct numerical simulations of the full Leggett-Takagi equation of motion for 
magnetization in 3 He-B and experimental data for magnetic field dependence of the catastrophic relaxation, 
that provide solid support of the theory of this phenomenon, presented in our 2006 JETP Letter. 
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Introduction. This Letter is inspired by critical com- 
ment on our JETP Letter "Solution of the problem of 
catastrophic relaxation of homogeneous spin precession 
in superfluid 3 He-B" The self sustained and long- 
lived spin precession with the coherent phase across the 
whole precessing domain - the Homogeneously Process- 
ing Domain (HPD) - is a unique feature of superfluid 
3 Hc-B. This phenomenon bears all the ingredients of 
spin superfluidity py. Later it was found that at low 
temperatures, below ~ (0.4 — 0.5)T C , the spin super- 
fluidity experiences abrupt instability, called the catas- 
trophic relaxation 0J. 

After many attempts it was finally recognized that 
the origin of the catastrophe is the parametric (Suhl) 
instability [SJ [2] ■ Two competing contributions to the 
increment of the parametric instability were suggested: 



V(r)=V w (r) + V BW {r) . 



(1) 



The "intrinsic" contribution V^ F , suggested in the 
Sourovtsev- Fomin (SF) Ref. [5], is due to the anisotropy 
of spin- wave velocity and it comes from the region where 
the configuration of spin and orbital vectors S and L is 
canonical, the so-called the Brinkman-Smith (BS) mode 
of precession. Suggested in our Bunkov-L'vov-Volovik 
(BLV) Ref. [2] contribution V BIjV is due the oscillat- 
ing spin-orbit energy and it only comes from the region 
where precession deviates from the BS mode. 
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Comment PQ clarifies aspects of the problem that 
may cause misunderstanding and thus require more de- 
tailed explanation and even further development. This 
is the subject of our Letter, which is organized as fol- 
lows. Sections 1-3 are devoted to questions, arose in P, 
which arc mainly related to the problem of parametric 
instability in the case of the spatially inhomogeneous 
precession. In Sec. 4 we consider extension of our the- 
ory to other precessing states. In Sees. 5 and 6 we dis- 
cuss results of direct numerical simulations of the full 
Leggett-Takagi equations for magnetization and of lab- 
oratory study of the magnetic field dependence of the 
catastrophic relaxation, that provide additional support 
of our theory of this phenomena, presented in Ref. |2] , 
and demonstrates that the SF mechanism of the insta- 
bility 5 is in contradiction with the experiments. 
1. Stability of the basic reference state. This is 
the first question of Com. PP that is useful to clarify. 
As stated in pp, the spin-orbit potential energy as a 
function of S z has concave or the convex form shape, 
depending on orientations of L and S and therefore in 
some cases the spatially homogeneous precession can 
be unstable 0. For example, the conventional homo- 
geneous precession is unstable in 3 He-A [7j. However, 
this is not applicable to a "texture" with some profile 
of the orientations of the orbital and spin moments, L 
and S, that realizes the minimum of the sum of the 
gradient and the spin-orbit energies (see Ref. (SJ for the 
texture in slab geometry and Ref. [H] for the texture 
under discussion). The texture is stable with respect 
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to perturbations with the same rotation frequency l>l- 
This texture, shown in Fig. 1 of our Ref serves as the 
basic reference state, stability of which with respect to 
parametric excitation of the spin waves with frequency 
uj s (±k) = ujl/2, was studied in our Ref. 0. Notice 
again, that the texture is stable with respect to perturba- 
tions with the frequency ujl- This is the answer to the 
first question of pQ : "Why the same mechanism which is 
responsible for the instability of the uniform precession 
in 3 He-A can be disregarded?" . 

2. Parametric instability in spatially inhomoge- 
neous media. The second statement of Com. PQ is 
trivial: because the spatial variations of the paramet- 
ric instability increment V "is not small ... , two ways 
of averaging [ exp((V)t) and (exp(Ui)) ] can give very 
different results". The relevant question here: "What 
is the adequate way of averaging of the instability in- 
crement (if it exists)?" To answer this question one has 
to describe parametric instability in spatially inhomoge- 
neous media. 

This problem has been discussed in various physical 
situations and is presented in many books, for example 
in Sec. 6.5.2 of monograph [TJ]|, where explicit expres- 
sions for the instability increment for different types of 
spatial inhomogeneities are derived. Without going into 
details we can say that the main physical message of 
this study is that for weakly decaying parametric waves 
(which is our case) the threshold of the parametric insta- 
bility can be estimated with a good accuracy from the 
total energy balance in the sample (cell, in our case). 
Namely, at the threshold the total energy influx into 
the system 

W+(x J V[r,k(r)]dr, (2) 

has to be equal to the total rate of the energy dissipation 
W- . With the same prefactor as in Eq. @ 

W„ oc J j[k(r)]dr, (3) 

where "f(r) is the damping rate of the parametric waves. 

The physical reason is that under wave propagation 
in weakly inhomogeneous media its frequency serves as 
the adiabatic invariant. This means that the wave fre- 
quency is independent of the position, while the wave- 
vector k(r) changes accordingly to the dispersion law 
ui(k,r) = const. Because of that, even propagating in 
inhomogeneous media, the waves with 

w[fc(r),r]=Wi/2, (4) 

stay in parametric resonance with the pumping (homo- 
geneous in frequency, but spatially inhomogeneous spin 



precession). In other words, if the mean- free pass of the 
wave exceeds the cell size, the entire cell can be con- 
sidered as the resonator with some parametrical mode, 
which is locally close to the planar spin wave. Under this 
conditions, the threshold can be well estimated from the 
condition: 

J V[r,k(r)}dr = J j[k(r)]dr . (5) 

Accuracy of this estimate is about ±(10-=-30)%, and 
related with ignoring in Eq. J^J) the spatial variation 
of the actual profile of the parametric waves (for more 
details, see i.e. Ref. \Tty). 

In experimental conditions, (see Fig. 1 in our 
Ref. 0) one roughly says, that the fraction of the cell 
volume (near the wall), with essential deviation from the 
planar geometry is about of a half of the total volume. 
This is the volume, where the energy pumping is domi- 
nating [i.e. the integral in the LHS of Eq. (JSJ]. In this 
region, in agreement with the the third comment of PP, 
the spin-wave vectors are relatively small and thus one 
can neglect the wave damping. The wave damping [the 
integral in the RHS of Eq. JSJ] is dominated by the cen- 
tral part of the cell, where wave vectors are sufficiently 
large. Accordingly, the total energy balance in the cell 
is kept by the spatial energy flux from the near-wall 
region to the central part of the cell. Under this condi- 
tions, the threshold Eq. JSJ can be roughly rewritten in 
the simple manner: 

Vv maxV[r, k(r)] ~ V 7 max7[fc(r)] , (6) 

in spite of the fact, that the energy pumping and the 
energy damping are dominating in the different regions 
of the cell. Here V 7 and Vv are the corresponding effec- 
tive volumes which depend on texture. As follows from 
the analytical and numerical simulations, these volumes 
do not depend on the value of magnetic field (see be- 
low) and in the considered cell are roughly equal. This 
means that 

maxU[r, fc(r)] ~ max7[fc(r)] , (7) 
estimate, used in our Ref. 0. 

3. Wave damping in the near-wall region. In the 

third complain of pQ we have been instructed how to 
estimate from Eq. J3J wave vector of k(r). Presented 
estimate shows that near the wall wave vectors and, 
consequently, wave damping is small with respect of the 
central region with the conclusion that we miss a factor 
about 20 in our estimate of the parametric threshold. 
The misunderstanding of Ref. pQ in this point is related 
to the question: ''How the value of the wave damping 
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in the near-wall region effects on the threshold of the 
parametric excitation in the cell?" We hope that our 
explanation in previous Sec. 2 is clear enough to re- 
alize that the energy balance in the system of weakly 
decaying waves has to be discussed rather globally, for 
the entire cell, then locally, point-wise, as presented in 
Comment pQ . Therefore, the estimates of the spin- wave 
vectors in the near- wall region, made in Com. pQ, be- 
ing reasonable itself, are irrelevant to the problem at 
hand. The same we can say about the statement, made 
in , concerning "additional factor w 20 in the RHS of 
Eq. (28)" in our Ref. As we explained in Sec. 2, the 
possible inaccuracy of our estimate of the threshold 
does not exceed (20 -j- 30)%. 

4. Generalization on other precessing states. To 

conclude the subject of the spatial inhomogeneity we 
mention, that the precession under discussion - the 
HPD - is a very specific precessing state due to its 
unique symmetry. Only in the case of the Brinkman- 
Smith mode the spin-orbit interaction does not con- 
tribute to the amplitude of the parametric instability V , 
and thus the parametric excitation by the BLV mecha- 
nism requires the existence of the texture. 

However, there are many other modes of spin preces- 
sion in 3 He, for which the spin-orbit potential energy as 
a function of S z has also the concave form shape. These 
are: the so-called HPD2 in 3 He-B JT]; the special mode 
of precession in 3 He-A ^2]; the precession at one half of 
equilibrium magnetization and at almost zero magneti- 
zation observed in 3 He-B ^2]; etc. As distinct from the 
HPD based on the Brinkman-Smith mode, in all these 
modes the spin-orbit energy is oscillating, and thus it 
produces the non-zero contribution to V even in case 
of the spatially homogeneous precession. This means 
that for all these modes of spatially homogeneous pre- 
cession, the BLV mechanism, based on the spin-orbit 
energy, will compete with the SF mechanism even at 
moderate magnetic fields. 

Notice, that some of the predicted HPD modes still 
have not been observed or identified. Now with our 
knowledge of the parametric instability mechanisms 
and their dependence on different parameters, we are 
able to find the region of parameters (magnetic field, 
temperature, angles of precession, superfluid velocity, 
etc.) where all the conditions for the stability of 
precession are satisfied. 

Now we are coming to the philosophical discussion, 
open in Com. pQ of "validity of an idealization for a 
particular experimental set-up" and of acceptance of 
"generalization of theory in a way, that makes possi- 
ble its application to wider class of experimental condi- 



tions" . Obviously, a success of "idealization" or "gener- 
alization", used in a theory, depends on an experience, 
physical intuition, taste and courage of its authors, that 
are different for different investigators. We hope that we 
made it clear in Sees. 1-3, that what it seems for the 
author of Com. pQ "ambiguous assumptions and non 
justified approximations" have indeed clear sense and 
well established physical background. 

Nevertheless we are happy, using the discussion with 
Com. pQ, to provide our theory of catastrophic relax- 
ation [5] with additional support from the direct nu- 
merical simulations of full Leggett-Takagi equations [H] 
(in the next Sec. 5) and especially from the laboratory 
study of magnetic field dependence of the temperature 
of catastrophic relaxation in 3 He, Sec. 6. Note that 
most of the complications related to the spatially inho- 
mogeneous precession disappear at high magnetic field, 
where the wave vector k is almost constant across the 
cell. 

5. Discussion of numerical experiment. One- 
dimensional numerical simulations of the Leggett- 
Takagi equations were made in Ref. 14 . It was found 
that in the viscous limit, i.e. at sufficiently high 
temperature, the Brinkman-Smith configuration takes 
place in the central part of the cell being disturbed 
by the periodic perturbations in angle 9 coming from 
the peripheral region, where the vector L considerably 
deviates from the Brinkman-Smith configuration. This 
observation confirms our statement, made in Sec. 1, 
that the spatial inhomogeneous texture is stable, pro- 
viding minimum of the spin-orbit and gradient energies 
at given boundary conditions and can not be destroyed 
by any relaxation processes. 

At further cooling of the cell, the BS configuration is 
more and more disturbed, and at some at some thresh- 
old temperature T t h the catastrophe occurs: the expo- 
nential growth of standing spin waves [with frequency 
lul/2] starts in the peripheral region and propagates to 
the central part, taking energy from the homogeneous 
precession and finally destroying it. This picture is in 
full agreement with the parametric instability of the in- 
homogeneous precession, studied in Refs. |2J|S]- 

Let us remind that there are two competing contri- 
butions to the increment of the parametric instability 
in Eq. JJ). The "intrinsic" contribution V SF comes from 
the whole sample; and the contribution V BLV which is 
due the oscillating spin-orbit energy and it only comes 
from the regions where the configuration deviates from 
the Brinkman-Smith mode. In the one-dimensional ge- 
ometry used in the numerics |14) . the SF contribution is 
suppressed in the central region because the increment 
is zero for the waves radiated in the direction perpendic- 
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The experimental data for the temperature of catas- 
trophic relaxation as a function of the NMR frequency 
at 31 bar ( o ), left scale, and the estimated value of the 
inverse spin diffusion coefficient, ( ) right scale. Solid 
line - BLV theoretical estimation |5J, dashed line - SF 
results 0. 



ular to L (cf. Eq.(26) of Ref. and Eq.(25) in Ref. 0. 
This means that while these particular numerical sim- 
ulations support the process discussed by the BLV the- 
ory 0)j they cannot either support or disregard the 
relevance of SF contribution at given parameters of the 
simulations. 

6. Magnetic field dependence: SF vs. BLV 
competition. To find out which contribution is domi- 
nating in real experimental situations, we will compare 
both contributions to the increment JIJ [21 directly 
with laboratory experiments. The BLV arguments |2] 
demonstrate that in the moderate magnetic fields the 
V BLV contribution is dominating. Both groups P El 
agree, that according to the calculations of Ref. [5] , the 
SF contribution V SF dominates at large magnetic field, 
when ojl 3> ^l- To see that this is really the case, 
we examine here experimental data on magnetic field 
dependence, that include the region of higher magnetic 
fields. 

The most relevant for this goal data are provided by 
the Cornell ^H] and Grenoble ^B] experiments, where 
the temperature of catastrophic relaxation has been 
measured as a function of magnetic field, see Fig. 1. 
Open circles in Fig. 1 demonstrate the temperature 
T ca .t{^L) of the catastrophic relaxation as a function 
of the magnetic field (or Larmor frequency u>l) at 31 
bar ^5J. The main effect of the temperature is to pro- 
vide the dissipation which damps the parametric in- 
stability via the spin-diffusion mechanism [5J. That is 



why the relevant physical quantity is the spin-diffusion 
coefficient D(T,ujl), rather than the temperature it- 
self, and we must convert the temperature to diffu- 
sion. According to Ref. ^7], in the ballistic regime 
considered here when u>lt(T) 3> 1, the spin diffusion 
D(T,uil) ~ l/wx,r(T), where r(T) is a quasiparticles 
scattering time, which grows exponentially with cool- 
ing. The experimental temperature dependence of the 
spin diffusion is obtained from measurements of diffu- 
sion through the HPD boundary 18 . 

Combining these results with T ca t(u>L,) from Fig. 1 
one obtains the spin diffusion at which the catastrophic 
relaxation occurs as a function of w^: solid circles in 
Fig. 1 show D~Jj.(uiL). The experimental points fit ex- 
cellently the cubic dependence 
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D. 



(8) 



obtained within the BLV mechanism {see Eqs. (23) and 
(24) in Ref. g]}. 

For the quantitative comparison of the experiment 
for ujl dependence of l/D cat with Eq. JSJ we shall 
take the experimental values of the parameters f^, cm 
and use the parameter a, which characterizes the L- 
texture, as a fitting parameter. In the considered re- 
gion of temperature all these parameters are slow func- 
tions of temperature. We shall use the data at 1 MHz 
and T = 0.35T C , where reliable experimental data exist. 
From Ref. we can estimate fl 2 L = 10 11 Hz 2 ; while 
from measurements of different modes of HPD oscilla- 
tions |20j with pressure scaling by Fermi velocity, we 
find c| = 1.5 x 10 6 cm 2 /s 2 . By introducing these val- 
ues into Eq. JHJ, one obtains a = 0.07, which is in good 
agreement with the theoretical estimation of a w 0.1 for 
6 bar, made in Ref. 

In the current consideration it is important that the 
parameter a does not depend on magnetic field. The 
reason for that is that according to the Ref. [S], the 
characteristic length scale of the near wall region is 
~ c/ yjuj{ui — uil)- It is typically about c/CIl, and thus 
does not depend on u>l- 

Next step is to compare the experimental results 
with the SF contribution 151. It reads: 



d: 



(9) 



Here we accounted for that A max near the HPD bound- 
ary is equal 0.016o;l and therefore the field dependence 
of \/D cat should be linear. As one see in Fig. 1 this 
clearly contradicts to the experiment. Furthermore, if 
we plot the value of l/D cat , we find (Fig. 1) that the 
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SF result agrees with the experiment only in the re- 
gion of NMR of 500 - 1000 MHz. At higher fields the 
theoretical value of diffusion, at which catastrophic re- 
laxation with the SF mechanism should occur, definitely 
disagrees with the experiments. 

Discussion. According to our theoretical analysis, 
in moderate magnetic fields, i.e. at ujl smaller than 
about 1 MHz, the BLV contribution the parametric in- 
stability of HPD is dominating. At these fields the ef- 
fect of spatially inhomogeneity on the BLV mechanism 
is essential, and we clarified in Sees. 2 and 3 the corre- 
sponding points, rouse in Com. pQ. Most of the compli- 
cations related to the issue of spatial inhomogeneity do 
not arise at high magnetic fields, where lot, ^> Vlr,, and 
the wave vector k is (almost) homogeneous along the 
cell. However, we expected that at such a high field the 
BLV contribution is subleading , while the SF contribu- 
tion dominates, if lo\ > 10 tt L . The same opinion was 
expressed in Com. pQ, where it was stressed that the 
stronger magnetic fields have to be used for the experi- 
mental investigation of the "intrinsic" SF mechanism of 
catastrophic relaxation. 

On the contrary, the surprising experimental fact is 
that the magnetic field dependence of the catastrophic 
relaxation demonstrates that even up to a rather high 
field, when uj\ ~ 100 Q\, it is still quantitatively and 
qualitatively described by the BLV contribution to the 
parametric instability. Moreover, the magnetic field de- 
pendence is in striking disagreement with the SF con- 
tribution. The matter of the fact is that the SF contri- 
bution does not show up in these experiments at all. 

We do not think that explanation of this fact is re- 
lated to possible calculational mistakes in the analytics, 
presented in Ref. [HJ. We feel that reason(s) for the ob- 
vious qualitative disagreement between the SF analysis 
and experiment is deeper and may be related to possi- 
ble violations of the SF approach in the region of large 
k vectors. 

To make long story short, the BLV mechanism of the 
parametric instability gives quantitative description of 
the present experiments for any values of col used, while 
the SF contribution, that seems to be the leading one 
(for large ujl), is absent in the experiments by unknown 
reason. This is the main puzzle; its solution requires the 
further theoretical, numerical and experimental efforts, 
that are beyond the scope of this Letter. 

Acknowledgements. We acknowledge the critical 
comments by LA. Fomin [Q, that inspire us for further 
clarification and development of our theory of catas- 
trophic relaxation in superfluid 3 He-B, presented in this 
Letter. This work was done as the result of collabora- 
tion in the framework of the ESF Program COSLAB, 



the Large Scale Installation Program ULTI of the Euro- 
pean Union (contract number: RITA-CT-2003-505313) 
and the project ULTIMA of the "Agense National de 
Recherche", France. The work was also supported in 
part by the Russian Foundation for Fundamental Re- 
search and the US-Israel Binational Science Foundation. 



1. LA. Fomin, Comment on "solution of the problem of 
catastrophic relaxation of homogeneous spin precession 
in superfluid 3 He-B" , |cond^mat/0606760] 

2. Yu. M. Bunkov, V.S. L'vov and G.V. Volovik, Solu- 
tion of the problem of catastrophic relaxation of homo- 
geneous spin precession in superfluid 3 He-B", Pis'ma 
ZhETF, 83, 624-629 (2006 ) [JETP Lett. 83, 530-535 
(2006)], cond-mat/0605386 

3. See review: Yu.M. Bunkov, "Spin Supercurrent and 
Novel Properties of NMR in 3 He", in "Progress in 
Low Temp. Physics" 14, 69, ed. W. Halperin, Elsevier 
(1995). 

4. Yu.M. Bunkov, V.V. Dmitriev, Yu.M. Mukharsky, 
J. Nyeki and D. A. Sergatskov, Europhys. Lett., 8, 645 
(1989). 

5. E.V. Sourovtsev, LA. Fomin, Pis'ma ZhETF, 83, 479 
(2006). 

6. LA. Fomin, ZhETF, 78 2392, (1980) [JETP 51 1203 
(1980)]. 

7. LA. Fomin, Pis'ma ZhETF 30, 179 (1979) [JETP Lett. 
30, 164 (1979)]. 

8. L.R. Corruccini and D.D. Osheroff, Phys. Rev. B 17, 
126 (1978). 

9. Yu. M. Bunkov, O. D. Timofeevskaya and G. E. Volovik, 
Phys. Rev. Lett. 73, 1817 (1994). 

10. V.S. L'vov, Wave Turbulence Under Parametric Excita- 
tion, (Springer- Verlag, 1994) 

11. Yu. M. Bunkov and G.E. Volovik, ZhETF 103, 1619 
(1993) [JETP 76, 794, (1993)]. 

12. Yu. M. Bunkov, and G.E. Volovik, Europhys. Lett. 21, 
837-843 (1993). 

13. V.V. Dmitriev, I.V. Kosarev, M. Krusius, D.V. Ponarin, 
V.M.H. Ruutu, and G.E. Volovik, Phys. Rev. Lett. 78, 
86 (1997). 

14. Yu. M. Bunkov, V.L.Golo, O. D. Timofeevskaya Czech. 
J. Phys., 46, 213 (1996), J. Low. Temp. Phys., 137, 625 
(2004). 

15. DA. Geller, D.M. Lee, Phys. Rev. Lett. 85, 1032 (2000). 

16. Yu. M. Bunkov, J. Low. Temp. Phys. 135, 337 (2004). 

17. D. Einzel, J. Low. Temp. Phys. 84, 321 (1991). 

18. Yu.M. Bunkov, V.V. Dmitriev, A. M. Markelov, 
Yu. M. Mukharsky and D. Einzel, Phys. Rev. Lett. 65, 
867 (1990). 

19. P.J. Hakonen, M. Krusius, M.M. Salomaa, et al., J. Low. 
Temp. Phys. 76, 225 (1989). 

20. Yu.M. Bunkov, V.V. Dmitriev and Yu.M. Mukharsky 
Physica B 178, 196 (1992). 



